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We present measurements of the energy relaxation length scale i in two-dimensional electron gases 
(2DEGs). A temperature gradient is established in the 2DEG by means of a heating current, and 
then the elevated electron temperature Te is estimated by measuring the resultant thermovoltage 
signal across a pair of deferentially biased bar-gates. We adapt a model by Rojek and Konig 
[Phys. Rev. B 90, 115403 (2014)] to analyse the thermovoltage signal and as a result extract i, 
Te, and the power-law exponent a-, for inelastic scattering events in the 2DEG. We show that in 
high-mobility 2DEGs, i can attain macroscopic values of several hundred microns, but decreases 
rapidly as the carrier density n is decreased. Our work demonstrates a versatile low-temperature 
thermometry scheme, and the results provide important insights into heat transport mechanisms 
in low-dimensional systems and nanostructures. These insights will be vital for practical design 
considerations of future nanoelectronic circuits. 


There currently exist well-established methods to 
probe the low-temperature (low-T) electrical and ther¬ 
moelectric properties of two-dimensional electron gases 
(2DEGs). However, probing heat transport mechanisms 
in these systems has proven more challenging, primarily 
due to the lack of convenient low-T thermometers that 
couple directly to the electron gas. Conventional low- 
T thermometers such as germanium or ruthenium-oxide 
films are sensitive only to the lattice temperature Tl, 
the temperature of the crystal that hosts the 2DEG. At 
Tl < 1K the coupling between electrons and phonons 
becomes relatively weak, and therefore Te can differ sig- 
nihcantly from Tl. The electrical resistance of the 2DEG 
itself becomes insensitive to Te at these temperatures, 
since the majority of scattering events are from static 
impurities. These both therefore become ineffective at 
measuring Te in this regime. Accurate measures of Te can 
be obtained from the Coulomb-blockade characteristics of 
quantum dots [T] which are broadened at a finite Te. The 
weak localization characteristics of 2DEGs [5] can also be 
useful as they depend sensitively on the phase-coherence 
length which is Te-dependent. However, neither of these 
methods lend themselves easily to the measurement of 
spatial temperature gradients, which is required in order 
to measure the thermal conductivity k mil] or the energy 
relaxation length scale ^ in 2DEGs. 

Appleyard et al. m showed that the diffusive compo¬ 
nent of the thermopower S can be used to detect dif¬ 
ferences between Te and Tl. Here S = Vth/AT, where 
Vth is the thermovoltage developed in response to a tem¬ 
perature difference AT. To do this they measured Vth 
across a pair of quantum point contacts (QPCs) with a 
heated electron gas between them, as was done similarly 
in Ref. [5 ]j more recently in Ref. [7|. Then AT was 
estimated as Vth/S, where S was obtained using the Mott 
relation [8]- 
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EIG. 1. Device geometry, (a) False-colour SEM image of 
the device used in our experiments. The first bar-gate ther¬ 
mometer (BGT) is outlined by white dashed lines. The three 
contacts at the top of the device allow a four-terminal mea¬ 
surement of the heating element resistance, (b) Schematic 
representation of a BGT which consists of gated (II and III, 
yellow) and ungated (I and IV, blue) regions of 2DEG, termi¬ 
nated by an ohmic contact. The orange section between the 
gates illustrates a region of ‘hot’ electrons. 


Here fee is the Boltzmann constant, e is the electron 
charge, h = /i/27r with h being Planck’s constant, a is 
the electrical conductivity, E is the total energy and n is 
the chemical potential. This technique has been used to 
measure k in quantum wires mm and energy-loss rates 
in 2DEGs |10| . Chickering et al. E] showed Eq. to be 
broadly valid in gated regions of a 2DEG between 0.8 K 
and 2 K and therein suggested the possibility of using a 
symmetric pair of bar-gates as a low-T thermometer for 
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the electron gas. This method was recently employed to 
measure S in mesoscopic 2DEGs Usefully, the 

relatively large size of the bar-gates eliminates the need 
for electron beam lithography which simplifies the fabri¬ 
cation process. However, it was noted in Ref. m that 
the data systematically deviated from the Mott predic¬ 
tion. Rojek et al. na attributed these deviations to the 
spatial extent of the bar-gate thermometers (BGTs) be¬ 
ing comparable to the energy relaxation length (. in the 
2DEG and developed a model to account for this. In this 
work we adapt the model developed by Rojek and Konig 
to refine the analysis of the signal produced by a BGT, 
and to make an accurate measurement of 1. 

Figure [T^ shows a false-colour SEM image of a typical 
device. Please see the supplementary Material [TB] for 
wafer and fabrication details. The device consists of a 
heating element and a longitudinal strip of 2DEG (of 
width 100 pm and length Lstrip = 1 mm) which together 
form a ‘T’ shape, with three BGTs along the strip. The 
first BGT is at a distance Lti = 200 pm from the heating 
element. A top-gate sits over the heating element and the 
strip and is used to tune the electron density n in these 
regions. This design minimizes any power reflection at 
the interface between the heating element and the strip. 
The strip is terminated by a large ohmic contact. 

Throughout the experiment all ohmic contacts through 
which no current is passed were assumed to be at Tl, 
since they are in direct thermal contact to the mixing 
chamber via the measurement wiring. Figure shows 
a schematic of a single BGT. It is a symmetric structure 
consisting of two arms flanking the 2DEG strip. The 
left (right) arm is formed by a gated region labelled II 
(III), followed by an ungated region labelled I (IV), and 
terminated by an ohmic contact. The lengths of the gated 
and ungated regions of the arms are Lg = 150 pm and 
Lug = 455 pm, respectively. 

The experiment involves passing a current Jh at fre¬ 
quency / = 10 Hz through the heating element which 
Joule heats the electron gas with a power of I^Rh., where 
i?h is the four-terminal electrical resistance of the heating 
element. This establishes a temperature gradient along 
the length of the 2DEG strip. The thermovoltage Uth 
generated across a thermometer in response to Jh is de¬ 
tected at 2/ using a lock-in amplifier. Figure shows 
4th = ha — Vb (see Fig. [^) from the first BGT while the 
gate voltage (I 4 ) on gate H is swept. An almost identical 
result is obtained when gate HI is swept, except that Uth 
is opposite in sign. The thermovoltages developed across 
the second and third BGTs were found to be negligible 
for even the largest used heating currents. Each section 
of the BGT from I to IV, contributes to Vth- 

Vth = ^ SpATp. (2) 

p 

Here p goes from I to IV, denoting each section shown 
in Fig. [^, ATp is the temperature drop across section p, 
and Sp for each section is given by the Mott formula for 


h(V) 
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FIG. 2. The fit shown in this figure is based on a simple 
model where the hot electrons are assumed to relax to Tl 
within the gated region. This model ignores any effects of the 
temperature profile along the BGT, and therefore the only 
fitting parameter is Te. The fit is clearly inadequate with the 
fit systematically deviating from the experimental data. This 
kind of deviation was found to be typical across the parameter 
space of Jh and Kg explored in this experiment. 

a non-interacting 2DEG: 

^ _ TTkB^Tpm (1 -b ae) 

^ ~ 3eJ2 Up ' ^ ^ 

Equation!^ is obtained by substituting the Drude ex¬ 
pression a = npe^Tg/m for the electrical conductivity into 
Eq. Here rip is the 2D number density of charge car¬ 
riers in section p, is the Drude elastic scattering time, 
m is the effective mass of the charge carriers (= 0.067me 
in GaAs-based 2DEGs, with nie being the bare electron 
mass), ae = inp/Te){dTe/dnp), and Tp is the average elec¬ 
tron temperature in region p. 

Figure also shows the best fit to the data by assum¬ 
ing that hot electrons relax to Tl within a distance Lg 
and therefore over the gated region, regardless of the car¬ 
rier density rig beneath the gate. In this simple picture 
the temperature profile along the BGT arm is not rele¬ 
vant, and the only unknown is between the BGT arms. 
However the quality of the fit is clearly inadequate and 
this is found to be the case across the parameter space 
of Jh and Ufcg explored in this experiment. 

We first describe why it is essential to consider the con¬ 
tribution from all the regions p = I to IV towards Kh- 
At low-T and especially in high-mobility 2DEGs, the en¬ 
ergy relaxation length i over which hot electrons relax 
through inelastic processes can greatly exceed the mean- 
free path of electrons [Him. The dependence of I on the 
inelastic scattering time r, is given by £ = where 

D = is the diffusion constant of the electrons, 

and Up is the Fermi velocity, r, has a power-law depen¬ 
dence on n: r, = Ti_o(^/''^o)“S where the subscript 0 de¬ 
notes the values in the ungated 2DEG. Thus, the distance 







3 



FIG. 3. Fits to the data using Eqs.[^|^and[^ The graphs show Vth as a function of the swept gate voltage 14 Ug, the 
corresponding carrier density under the swept gate. The different traces in (a) and (b) represent different heating currents 7h 
and different top-gate voltages Vtg respectively. The fits are markedly improved over, for example, that shown in Fig. 
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FIG. 4. (a) The elevated electron temperature ATe = Te — TL extracted from the model is seen to increase as a power-law of 7h. 
The broken line shows an 7^ trend, which would be the case if ATf. increased linearly as a function of applied heating power. 
The sub-linear dependence on power indicates that there are heat losses to the surrounding lattice, (b) As Vtg becomes more 
negative, Te at the first thermometer decreases rapidly, indicating that the electrons thermalize over a much shorter distance. 


over which electrons lose their excess energy in the BGT 
arm is crucially dependent on Up. The ohmic contact en¬ 
forces Te = Tl at the 2DEG to ohmic contact interface, 
and this needs to be taken into consideration if i is com¬ 
parable to Tg-|-Tug- Therefore in all the above situations, 
both the gated and ungated regions contribute to Vth in 
a manner dependent on Tg at the centre of the BGT, Tg 


at the interface of the gated and ungated regions, £, and 
tti, all of which need to be estimated self-consistently. To 
do this we adapt the model used in Ref. [TS] (described 
in the Supplementary Material [E]). The resulting ex¬ 
pression for ATe = Te — Tl at the junction of the gated 
and ungated regions ATe(Tg), reads: 


ATe(Tg) = ATe 


_zsinh(Tug/7o)_ 

cosh(Tug/7'o) sinh(Tg/7) -|- z sinh(Tug/^o) cosh(Tg/7) 


(4) 


Here £o is the energy relaxation length at n = no, and 
z = Within the framework of Ro- 

jek et al. ’s linearized model US], Tp = Tl- Equation]^ 


provides an expression for ATp, which when substituted 
together with Eq. in Eq. results in an expression 
for Vth as a function of no, n in each gated region, £o, 
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Qfe, cti and Te- We measure ng and n(Vg) in the device 
by observing 4g-dependent edge-state reflections in the 
quantum Hall regime |16j . and is extracted from the 
dependence of tr on n and turns out to be ~ 0.89 over the 
relevant range of n. This leaves three unknowns, namely 
£o, Te, and a,, which are used as fitting parameters. Im¬ 
portantly though, by fitting to several complementary 
data sets whilst varying different experimental param¬ 
eters such as Ih and Vtg, we are able to considerably 
reduce the uncertainty in these three fitting parameters 
(see Supplementary Material m)- 

Figure shows Vth against rig for varying Ig which 
will produce different ATg. Figure shows Vth against 
rig but for varying Vtg which will vary £ in the 2DEG strip 
as well as ATe via its effect on Rh- Clearly, the model 
produces excellent fits to the data with no discernible 
systematic deviations. Similar data and quality of fit is 
obtained when the adjacent BGT is swept. 

The results of the fitting are the ATe for each R and 
Vtg in the two data sets. Figure]^ shows ATe as a func¬ 
tion of Ih on a log-log scale, and we find that ATe oc Ih^^- 
This sub-squared dependence is presumably due to a frac¬ 
tion of the power being lost to the lattice and to the ohmic 
contacts. The values of a-, and tg are found to be « 3.76 
and « 280pm, respectively. Figure]^ shows ATe as a 
function of Vtg after ATe has been scaled for the changing 
value of Rh- This has been done by applying a corrective 
factor of Rh,o/Rh{Vtg) so that any change in ATe is now 
due solely to a change in t. Here Rh{Vtg) is the electri¬ 
cal resistance of the heating element as a function of Vtg. 
Therefore the graph suggests that for Vtg < —0.11 V , the 
hot electrons completely relax within a distance Tti- 

The data therefore suggests that hot electrons ther- 
malize over macroscopic length scales (« 300 pm) in the 
ungated 2DEG, but that this length scale rapidly de¬ 
creases with n. This strongly justifies the need to ac¬ 
count for the ungated 2DEG arms when using BGTs at 
high Tig. This is also consistent with the negligible ther¬ 
movoltage detected across the second and third BGTs 
which are at distances of 500 pm and 800 pm from the 
heating element, respectively. While rig can certainly be 
lowered untill £ < Lg such that ATe(Tg) = 0, care must 
be taken to ensure that Eq. [^remains valid. Indeed, we 
have observed that when the 2DEG approaches the lo¬ 
calized regime (when a becomes < jli) the model is 
unable to fit the data satisfactorily. In this experiment 
n was conservatively limited to > 0.7 x 10^® m“^ corre¬ 
sponding to k-pl > 4 and Ts < 2, where fcp is the Fermi 
wave vector, I is the elastic mean free path, and rg is the 
interaction parameter defined as the ratio of the Coulomb 
energy and kinetic energy of the 2DEG. 

Importantly, the model only provides the value of £g 
and tti from which £{n) can be reconstructed. However, 
as argued in the previous paragraphs, this information 
is independently contained in Fig. [^. The Te-profile in 
the 2DEG strip is given by (see Supplementary Mate¬ 
rial [TB]1: 
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FIG. 5. Energy relaxation length i as a, function of carrier 
density n. The data in Fig. |4jo produces an independent mea¬ 
sure of i{n) which we find to be in excellent agreement with 
the model exponents. 


ATe{y) = ATe,o 


sinh [(Tgtj-ip y)/f*] 
sinh (Tgtj-jp /£^ 


(5) 


Here y is the spatial coordinate along the 2DEG strip 
(see Fig. Ef) and ATe^o = ATe{y = 0), the temperature 
elevation in the heating channel. Since the heating chan¬ 
nel is effectively at a constant temperature between pairs 
of points in Fig. S’: the ratio ATe(ni)/ATe(n 2 ) taken at 
y = Tti provides an implicit relation between £{ni) and 
£{n 2 ), where ni and n 2 are the respective carrier densi¬ 
ties. Therefore this allows the inference of £ from Te, and 
reconstruction of ^ as a function of n as shown in Fig. 
This is the main result of this study showing the depen¬ 
dence of the energy relaxation length scale on the carrier 
density. We stress that this is a direct measurement of 
£{n) using the pre-calibrated BGT, which we find to be in 
striking agreement with indirectly obtained dependence 
£{n) = £g{n/ng)^^~^°‘’‘~^°''R‘^, using the values of £g and a, 
derived from the model and the fitting. 

To summarize, we have demonstrated that BGTs are a 
versatile tool with which to detect elevated electron tem¬ 
peratures. As remarked earlier, BGTs are an attractive 
alternative to QPGs. QPG fabrication requires electron 
beam lithography and they come with the associated dif¬ 
ficulties of sub-micrometre devices such as sensitivity to 
electrical shock and vulnerability to disorder. In con¬ 
trast, BGTs can be macroscopic and therefore free from 
the above mentioned difficulties. 

While our manuscript seems to suggest that the trade¬ 
off between QPGs and BGTs is that detailed modelling 
is required to extract Te, it is important to note that 
modifying the device design such that Lg > £ significantly 
reduces the complexity of the analysis. On the other 
hand, the advantages of the employed device design are 
that it allows for the determination of £{n) using a single 
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BGT. We note that using a second BGT further displaced 
along the length of the 2DEG strip should also allow for 
such a measurement, so long as £{n) < Lt 2 - 
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Determining energy relaxation length scales in two-dimensional electron gases 

Supplementary Material 


I. WAFER, DEVICE FABRICATION, AND 
MEASUREMENT SYSTEM DETAILS 


The 2DEGs used in this study have a mobility of 
172m^/Vs at a carrier density of no = 1.45 x 10^^ m“^. 
For a top-gate voltage of Ug = 0 V, i?h = 230.9 ft which 
corresponds to a resistivity of 44.4 U/Sq. The typical 
ohmic contact resistance is « 200 U as determined by 
the difference between two and four-terminal resistances 
(corrected for known cryostat wiring resistances). There¬ 
fore the ohmic contacts should not effect thermalization. 
A wet etch is used to define the conductive mesa, and 
Au-Ge-Ni ohmic contacts and Ti-Au top-gates are then 
deposited by thermal evaporation. Photolithography is 
used at each stage to define the layer pattern in photore¬ 
sist on the substrate. Measurements on the devices were 
performed in a dilution refrigerator with a base temper¬ 
ature of 170 mK. A ruthenium oxide thermometer at¬ 
tached to the mixing chamber of the dilution refrigerator 
was used to measure Tl. 


II. THERMOPOWER DIFFUSION MODEL 

Rojek et al. consider a fixed electrical power input into 
a two-dimensional electron gas (2DEG) that results in an 
elevated electron temperature Tg = ATe+TL with respect 
to the lattice temperature Tl. In the steady-state, the 
power input to the 2DEG is balanced by the heat loss 
to the lattice. Then, one can derive from the continuity 
equation for heat and charge flow a diffusion equation for 
ATg{x), 


f' d 
AT - — — 

K OX 


dATg 
' dx 


= 0 . 


(SI) 


Here k is the thermal conductivity of the 2DEG, £ the 
energy relaxation length, and x is the spatial coordinate 
along which the temperature varies (for the temperature 
variation within the BGT, relevant for Eq. 4, this is the 
cc-direction while for the temperature variation within the 
2DEG strip, relevant for Eq. 5, this is the y-direction, see 
Fig. 1 of the main text). 

An applied gate voltage changes the electron density 
n which, in turn, modifies the transport properties, in 
particular the elastic and inelastic scattering times Tg and 
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Fig. S 1. Fits to a set of data showing the sensitivity of the fit to altered fitting parameters. 


Ti, the thermal conductivity k, and the energy relaxation 
length Their n-dependences are described by power 
laws with the following exponents: 

dlnTg/dlnn = a^, 
dlnTj/dlnn = Oi, 
dlnK/dlnn = 1 + tte, 
dln£/dlnn = (1 + Og + cxi)f2. 

For each region with constant electron density (the 
gated region II(III) and the ungated region I(IV) of the 
BGT, relevant for Eq. 4, as well as the 2DEG strip, rel¬ 
evant for Eq. 5), the general solution of Eq. [^is 

AT^{x)=ae-^/^+ be^/\ (S2) 

where the coefficients a and b have to be determined by 
boundary conditions. 


A. Derivation of Eq. 4 

For the BGT, the spatial position x is defined relative 
to the interface between the 2DEG strip and the gated 


region IL There are two regions of constant electron den¬ 
sity (the gated and the ungated one, I and II) which 
requires four boundary conditions to fully determine the 
temperature profile. One is given by ATe(O) at x = 0, 
the quantity to be measured by the BGT. Furthermore, 
the Ohmic contact at x = Lg-|-Lug fixes the electron tem¬ 
perature to the lattice temperature, ATe(Lg -|- Lug) = 0. 
Finally, at x = Lg (the interface between the gated (I) 
and the ungated (II) region), both the temperature ATg 
and the heat current —K{dATe/dx) have to be continu¬ 
ous. It is, then, a matter of solving a set of four linear 
equations to express the coefficients of Eq. in terms 
of Arg(O). However, the only relevant feature of the full 
temperature profile that enters the thermovoltage is the 
temperature at the interface between region I and II, 
ATg(Lg). Plugging x = Lg into the full solution im¬ 
mediately leads to the compact formula Eq. 4, where 
z = {K/i)/{Ko/^o) = (u/no)^^^“'’““‘^/^ originates from 
the boundary condition of the heat current to be contin¬ 
uous. 

In the limit of a large gated region, Lg £, the coeffi¬ 
cient b of the exponentially increasing term in Eq. [^van¬ 
ishes, ATe(Lg) goes to zero, and the ungated region does 
not influence the measure thermovoltage. For Lg < £, 
however, the ungated region becomes important. Ignor¬ 
ing it leads to an inadequate fit to the data, as shown in 
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Fig. 2 of the main text. 


B. Derivation of Eq. 5 

The modelling of the temperature profile along the y- 
direction in the 2DEG strip is analogous to that of the 
BGT. Since the electron density is constant in the 2DEG 
strip, only two boundary conditions are needed. They 
are given by fixing ATg to ATg^ at y = 0 and to 0 at 
the Ohmic contact y = Lstrip- Formally, determining 
the temperature profile along in the 2DEG strip is just 
a special case of the calculation for the BGT. Therefore 
Eq. 5 follows from Eq. 4 by performing the replacements 
z — > 1, fo f, Lg —>■ y, Tug —)■ Tstrip ~ y, and ATg — >■ 

^Tg^o. 


III. SENSITIVITY OF FIT ON PARAMETERS 

Figure S 1 shows the sensitivity of a model fit to the 
fitting parameters. The particular data used to demon¬ 
strate this corresponds to Ih = 600 nA at Vtg = OV, 
although similar results are obtained for all other values 
of Ih and Vtg used in this experiment. The optimum fit 
is contrasted to the resulting fit when a specific fitting 
parameter is changed by a given percentage. In panel 

(a) ATg was varied by 10% in each direction, in panel 

(b) and (c), Oi and £ were varied by 20% in each direc¬ 
tion respectively. The result is a clearly deviated and 
inadequate fit in each case. However these percentage 
changes are chosen such that the quality of the fit vis¬ 
ibly deteriorates. In fact the residual of the fit is more 
sensitive to changes in a single fitting parameter, for ex¬ 
ample, doubling for a 5% change in ATg. This highlights 
the technique’s usefulness as an electron thermometer. 



